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We study the large-scale inhomogeneity of the Universe based on the av- 
eraging procedure of Buchert and Ehlers. The generalized Dyer-Roeder equa- 
tion for the angular diameter distance of the inhomogeneous Universe is de- 
rived and solved for different cosmological models. We make a comparison 
of certain cosmic observables, such as the Hubble function, angular diameter 
distance, cosmological correction factor of homogeneous and inhomogeneous 
cosmological models, which are crucial ingredients in galaxy number counts 
and gravitational lenses. 

It seems that present cosmological data consisting of supernova searches, 
structure formation, gravitational lenses, etc. suggest a rather unexpected 
equation of state for the simple homogeneous Friedmann-Lemaztre (FL) Uni- 
verse with Q m ~ 0.3, Q\ ~ 0.7. However, one should be cautious even when 
performing analyses of the cosmological data (for example, for supernova 
data see the reanalysis in []!]). Further concern is the low statistical confi- 
dence of all extracted observables. Nevertheless, this phenomenological result 
initiated a lot of speculative work on the origin of the positive cosmological 
constant. 

Contrary to current investigations, there is only one derivation of the 
cosmological constant without any fine tuning or any obscure assumptions 
and this derivation is based on the Einstein-Cartan gravity Q. When the 
Universe is frozen at zero temperature T 7 = (spacelike infinity) and if all 
nonrelativistic fermionic matter is spinning matter, then owing to the same 



1 



coupling between spacetime curvature and mass density on the one side and 
spacetime torsion and spin density of matter on the other side, the following 
relations emerge: Q m = 2, = —10. It is important to strengthen 
that at T 7 = all spin densities of nonbaryonic and baryonic species act 
coherently to the total spin. The role of spin densities in the derivation of 
the primordial density fluctuation is described in Ref. ||. The existence 
of vorticity and acceleration 0, |3| or nonvanishing shear |4j] should be of 
fundamental importance for resolving certain cosmological problems, but the 
cosmological distance measures are not very much affected by relatively small 
deviations from the FL model . To conclude, let us notice that it is possible 
to construct a nonsingular and nonanomalous gauge theory with heavy and 
light neutrinos as cold and hot dark matter fermionic particles respectively, 
but without Higgs scalars and with no asymptotoic freedom in QCD |J. 

In this paper we want to reconcile current measurements with large-scale 
inhomogeneous cosmic models. We assume that the Universe is inhomoge- 
neous on large scales, thus dumpiness should not be characteristic only of 
small-scale structures, and that the Earth is placed in the region away from 
the centre of homogeneity. 

Starting with the historical work of Lemaitre, Tolman and Bondi (LTB) 
on inhomogeneous models, one can account enormous activity in this field . 
However, a recent attempt |3| to exploit LTB models, following perturbative 
calculations in Ref. ||, was not successful at large redshifts. 

The cellular structure of the Universe with a power-law distribution of 
matter of Ruflini et al flOfl , represents another approach to the problem of 
inhomogeneity. 



The work of Zalaletdinov |Tl[ is the first complete and consistent treat- 
ment of the averaged Einstein equations in an arbitrary Riemannian space- 
time. However, even if the vorticity and acceleration do not vanish, they are 
small deviations (which we neglect in this paper) from the FL geometry, thus 
we choose the following line element and spatially averaged scalar quantites 



of Buchert and Ehlers [12 



ds 2 = dt 2 - g^dX^X^ 

{m{t,x l )) v = ^- f d 3 xjv(t,x* 

V-n JT> 
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V v (t)= f Jd 3 X, J(t,X l ) = [det(g tJ )}^. 
Jv 

Introducing natural definitions 



a v (t) -- 



V 



{0)v = £ = A = is(t), 
V-d clv at 

M v = ( d^XJpit^X 1 ), 
Jv 

V V a V 

one can easily derive Raychaudhuri and constraint equations for averaged 



scalars from Einstein and conservation equations [12| [U| 



d v M- 



3— + ^G n ~P^y — A = Qx>, (1) 



3(^) 2 - SttG*-^- + W - A = -ig D , (2) 



Q 7y = l((e-(6) v ) 2 )- D -2(a 2 ) v , 



1Z = 1Z\, TZij = spatial part of the Ricci tensor. 



In addition, we need a condition of integrability of these equations [12 



(a b v Q v )- + aUa v (Jl} v y = 0. (3) 

Now we turn to the geometrical optics in order to derive the wanted 
observables of the inhomogeneous cosmic models. 
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Starting from the definition of the angular diameter distance and the 
focusing equation |T4}| : 



AttG n 
T 



_>| + _-(! + ,) VI = 0, (4) 



w = H ujqV, v = af fine parameter. 
We take into account that the shear of the congruence of photons vanishes 



for a spherically symmetric source 15 . Then, performing the Buchert-Ehlers 



spatial averaging and neglecting the spacetime shear in the averaged equation 
for the change of redshift with respect to the affine parameter: 

*(! + ,)»?-« (5) 
aw a>v{v) 

we arrive at the following equation for the averaged angular diameter 
distance: 



d 2 (D) v dz 2 d(D) v d 2 z AnG N 2 
+ ^z~d^ + ^lf {1 + Z) 

Assuming the smooth space dependence of mass density and angular di- 
ameter distance, one can factorize the last term of the preceding equation 

as 



( P D) V ~ p{t,X*)(D) v ~ ^(D) v , XeV, 

and we are left with the second-order generalized Dyer-Roeder equation 
with a suitable density normalization for the cosmological models: 



, ,d 2 D(z) . dD(z) , « ~ / \ 

[z) ^z^ + x{z) d^ + i^M*) = °' (6) 
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2 ; 



cf 2 2 . c?2;. 2 , 3 ^ 1 dhx>{z) . 
<iw; 2 c?w 1 + 2 a©(V) tte ' 
<idx> 1 



cfe #(z)(l + Z 

_ du(z) \ 1 



a v {z) l + z 



V(z) = l(l + z) 5 n m , A = 3Hfa A . 

Under adequate initial conditions [14]] (Adams-Bashforth integration method 
used) 



for integration from z\ to z 2 , z 2 > Z\ 

D(zi,z). \ 1 1 
D(zi, z\) = 0, V' > {z = z 1 ) 



dz H(zi) 1 + Zi ' 

we can calculate the averaged angular diameter distance for an arbitrary 
cosmological model Q m + Q A = 1- Needless to say, for the vanishing backre- 
action and curvature terms one recovers a homogeneous solution: 



H (l + z)J [Q m (i + C )3 + Q A ]i- 

Numerical evaluations and comparisons are performed in the following 
manner: (1) choose some functional form for the backreaction Q(z), (2) 
evaluate curvature term R(z) from the integrability condition, (3) choose two 
inhomogeneous models with Q m = 2, Q A = — 1 [Einstein-Cartan (EC) model] 
and fl m — 1 , Q\ — [Einstein-deSitter (EdS) model] and fit parameters 
of the backreaction and curvature to reach the input value for the Hubble 
constant Hq = a(z = 0)/a(z = 0) and to get the best fit of the D(z) of the 
inhomogeneous models to the D(z) of the homogeneous model with 
a m = 0.3 , tt A = 0.7 for z < 0.7. 
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One can easily check that for any > a > — 3 (we omit index V) 



Q(z) = Q a a (z), 
(K) v = R(t) = -a~ 2 (t) f ' dra- 4 (r)-^[a 6 (r)Q(r)i 

J to CLT 

R(z) = -^£[QW - ($r) 2 ^o)], a + -2, 

a + 2 a{z) 

R(z) = AQ(z) In a = -2, 

1 + z 

the curvature term power-law coefficient is within the same interval as 
that for the backreaction. Thus, this simple functional form garantees the 
homogeneity at very large scales | Q(z) |, | R(z) |<^C HQQ m ax>(z)~ 3 , z » 1 (it 
could also be worth including the Gaussian damping, if necessary). 

For definite numerical comparisons of the three cosmological models we 
use this Ansatz: 



Q(z) =Q [l + c a-\z)}, 

In order to match the input Hubble parameter at present time H , zq 
has to fulfil the following equation: 



Q(0) + R(0) = 
=>> 2(-l + 3c )(l + z ) 2 - 5c (l + Zq) + 3 = 0. 

Searching for the best fit values of c and Q , we find 

H = h u, u = lOOkms^Mpc" 1 , 
(EC) : ho = 0.6, Q = 3.789 w 2 , z = y/6/2 - 1, c = 0, 
(SriS) : h = 0.6, Q = 1-2632 w 2 , z = y/6/2 - 1, c = 0, 
(Horn) : /i = 0.6, Q m = 0.3, fi A = 0.7, Q = 0. 
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For z < 0.5, the angular diameter distances differ less than 5%. Figs. 
1 and 2 show the Hubble flow and the angular diameter distance for the 
three models, while Table 1 presents the cosmological correction terms for 



the difference in light travel time between gravitational lens images [16 



H At = Tf(6 1 

,Obs, V2,Obs, Zd, Z s 

D 2 (0,z d ) 



T = H n 



-il + Zd) 



z s - z d 



D 



D z s z d 
D(z d ,z s )D{0,z d ) 



D{0,z s 





z d = 0.5, z s = 1 


z d = 0.5, z s = 1.5 


z d =1, z s = 1.5 


z d = l, z s = 2 


T(EC) 


0.7550 


0.6686 


0.4977 


0.4226 


T(EdS) 


0.9559 


0.9350 


0.8775 


0.8471 


T(Hom) 


1.8543 


2.1098 


4.7948 


8.9155 



The galaxy number count depends essentially on the cosmological model 
through the comoving volume: 

dV(z) 



dz 



oc D z) 



We conclude with two observations: (1) For the large-scale inhomogene- 
ity of the Universe to be established one needs clear indications that none 
of homogeneous models can simultaneously fit and explain all data in cos- 
mography, structure evolution, gravitational lenses, etc. for all redshifts. At 
present it is difficult to make any conclusion what is the source of some re- 



cently found disparities in gravitational lenses JL7]] or galaxy evolution 
(2) If the inhomogeneity is established, then one could intend to model a 
backreaction function for the global fit of data, and ultimately attempt to 
derive it from the evolution of multicomponent imperfect fluid with photons, 
baryonic and nonbaryonic matter, with a necessary knowledge of all relevant 
masses and couplings. 



* * * 
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Figure 2: Solid, dashed and dot-dashed lines denote Hubble function H(z) 
in [u] for (EC), (EdS) and (Horn) models, respectively. 
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